Ab-initio study of the thermopower of biphenyl-based single-molecule junctions 
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Employing ab-initio electronic structure calculations combined with the non-equilibrium Green's 
function technique we study the dependence of the thermopower Q on the conformation in biphenyl- 
based single-molecule junctions. For the series of experimentally available biphenyl molecules, alkyl 
side chains allow us to gradually adjust the torsion angle tp between the two phenyl rings from 0° 
to 90° and to control in this way the degree of 7r-electron conjugation. Studying different anchoring 
groups and binding positions, our theory predicts that the absolute values of the thermopower 
decrease slightly towards larger torsion angles, following an a -I- 6cos^ ip dependence. The anchoring 
group determines the sign of Q and a, 6, simultaneously. Sulfur and amine groups give rise to 
Q,a,b > 0, while for cyano Q,a,b < 0. The different binding positions can lead to substantial 
variations of the thermopower mostly due to changes in the alignment of the frontier molecular 
orbital levels and the Fermi energy. We explain our ab-initio results in terms of a 7r-orbital tight- 
binding model and a minimal two-level model, which describes the pair of hybridizing frontier 
orbital states on the two phenyl rings. The variations of the thermopower with (p seem to be within 
experimental resolution. 



PACS numbers: 85.65.+h, 85.80.Fi, 73.63.Rt, 81.07.Pr 

I. INTRODUCTION 

Tailored nanostructures hold promise for improved ef- 
ficiencies of thermoelectric materials.^ For this reason 
there is a growing interest to gain a better understand- 
ing of the role of interfaces on thermoelectric properties 
at the atomic scale. Controlled metal-organic interfaces 
can be studied using single-molecule junctions, and re- 
cently the thermopower of these systems was determined 
in first experiments. 1^ While the thermopower (or Seebeck 
coefficient) of metallic atomic contacts was measured al- 
ready several years ago,-5- molecular junctions offer fas- 
cinating possibilities to adjust thermoelectric properties 
due to the control over chemical synthesis and interface 
structure. Ref. 2] and subsequent experimental studies 
thus explored the influence of different para meters on the 
thermopower, such as molecule lengthiSl^substituents,!^ 
anchoring groupSj^^HSl or electrode metaL^i^ 

On the theory side, the electronic contribution 
to the thermopower explains important experimental 
observations .'ii' We have shown recently that the ther- 
mopower of metallic atomic contacts, which serve as 
reference systems in molecular electronics, can be un- 
derstood by considering the electronic structure of dis- 
ordered junction geometries. Using molecular dynam- 
ics simulations of many junction stretching processes 
combined with tight-binding-based electronic structure 
and transport calculations, we found thermopower- 



conductance scatter plots similar to the low-temperature 
experiment.'^ Such a statistical analysis, although highly 
desirable for molecular junctions, is complicated by the 
time-consuming electronic structure calculations needed 
to describe these heteroatomic systems. Still, early stud- 
ies of the thermopower based on density functional theory 
(DFT) for selected geometries explained crucial trends, 
such as the dependence of the thermopower on molecule 
lengtlJ ^ or the influence of substituents and anchoring 
groups .'^^'^ Since the experiments on the thermopower 
of molecular junctions were all performed at room tem- 
perature until now, finite temperature effects may play 
a role. They can impact the thermopower by fluctu- 
ations of the j unction geometry and electron- vibration 
couplings .li^EIEni While the quantification constitutes an 
interesting challenge for future work, we will focus here 
on the purely electronic effects in static ground-state con- 
tact structures. 

An interesting aspect, not yet addressed in the ex- 
periments is the influence of conjugation on the ther- 
mopower Q. For the conductance, such studies were car- 
ried out by different groups with biphenyl molecules.l^^tl^ 
The torsion angle tp between the phenyl rings was 
adjusted stepwise by use of appropriate side groups. 
While such substituents may have a parasitic shift- 
ing effect on energies of current-carrying molecular or- 
bitals, the changes in conformation, which control the 
degree of 7r-electron conjugation, turned out to domi- 
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nate the behavior of the conductance!^ The systematic 
series of biphenyl molecules of Refs. fTSVETl uses alkyl 
chains of various lengths and methyl groups to adjust 
if and avoids strongly electron-donating and electron- 
withdrawing substituents. Hence it seems ideal for de- 
termining the influence of conjugation on thermopower. 

Theoretical work has considered the behavior of Q 
when if is changed continuously for the thiolat ed bip henyl 
molecule contacted to gold (Au) electrodes) ^^ * ^^ ! Both 
studies agree in the fact that Q is positive for all Lp. How- 
ever, while we predicted Q to decrease with increasing Lp 
based on DFT calculations and a 7r-orbital tight-binding 
model (TBM),J^work of Finch et aZ.l^^! suggested the op- 
posite for this idealized system. In this study we clar- 
ify this contradiction and demonstrate with the help of a 
two-level model (2LM) that for the off-resonant transport 
situation the absolute value of Q is expected to decrease 
when the molecule changes from planar to perpendicular 
ring orientation. This confirms our previous conclusions. 
More importantly, this work explores the possibility to 
measure the dependence of Q on Lp for the experimentally 
relevant family of molecules presented in Refs. [THH^T] 

Using DFT calculations of the electronic structure 
combined with the Landauer-Biittiker scattering formal- 
ism expressed with Green's function techniques, we de- 
termine the thermopower of biphenyl-derived molecules 
connected to gold electrodes. The molecules investigated 
are displayed in Fig. [l] Alkyl chains, one to four CH2 
units long, allow to change Lp gradually from 0° to 60°. 
To achieve ip « 90°, we included in addition M7, and as a 
reference also MO, the "standard" biphenyl molecule. For 
each of the molecules in Fig. [T| we will explore the three 
different anchoring groups sulfur (S), amine (NH2), and 
cyano (CN) in various binding geometries. 

This work is organized as follows. In Sec.|TT]we intro- 
duce the theoretical procedures used in this work. Sec. Ill 
presents the main results. We start by discussing mod- 
els to describe the Lp dependence of the thermopower, 
show the DFT-based results for Q, and provide further 
insights by discussing their relation to the predictions of 
the TBM and the 2LM. The paper ends with the conclu- 
sions in Sec. II VI 
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Figure 1: Chemical structure of the studied biphenyl 
molecules with X standing either for the S, NH2, or CN an- 
choring group. 



II. THEORETICAL METHODS 
A. Electronic structure and contact geometries 

We determine the electronic structure and contact 
geometries in the framework of DFT. All our calcula- 
tions are performed with the quantum chemistry package 
TURBOMOLE 6.3,1^ and we use the gradie nt-corrected 
BP86 exchange-correlation functional.l^^l^si -por the ba- 
sis set, we employ def2-SV(P) which is of split-valence 
quality with polarization functions on all non-hydrogen 
atoms.l^Sl For Au an effective core potential efficiently 
deals with the innermost 60 electrons,^^ while the ba- 
sis set provides an all-electron description for the rest of 
the atoms in this work. 

The contact geometries for the S-terminated molecules 
are those of Ref. |28j For the NH2 and CN anchors we pro- 
ceed as described in Refs. I19I28I and use for consistency 
the electrode geometry from Ref. [28] 



B. Charge transport 

We determine charge transport properties within the 
Landauer-Biittiker formalism. The transmission func- 
tion t{E), describing the energy-dependent transmission 
probability of electrons through the nanostructure, is 
calculated with non-equilibrium Green's function tech- 
niques. The Green's functions are constructed by use of 
the DFT electronic structure as obtained for the ground- 
state molecular junction geometries. A detailed descrip- 
tion of our quantum transport method is given in Ref. 1291 

The thermopower at the average temperature T is 
defined as the ratio of the induced voltage difference 
AV in the steady state and the applied temperature 
difference AT between the ends of a sample, Q{T) = 
— {AV/ AT)\j^Q. In the Landauer-Biittiker formalism 
the electronic contribution to the thermopower can be 
expressed adP^l 



QiT) 



Ki{T) 
eTKo{T) 



(1) 



with KniT) = / dET{E){E - fi)''[-df{E,T)/dE], the 
absolute value of the electron charge e = |e|, the Fermi 
function f{E,T) = {exp[{E- iJ,)/kBT] + l}-^ , the Boltz- 
mann constant ks, and the chemical potential fi^ Ep ^ 
— 5 eV, which approximately equals the Fermi energy Ep 
of the Au electrodes. At low temperatures, Eq. ([l]) sim- 
plifies tcP 



Q{T) = ~q{T) 

with the prefactor q{T) 
early on temperature. 



dET{E) 
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t{E) 
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Figure 2: (Color online) Analyzed types of junctions, shown for M3. For S anchors we consider hollow, bridge, and top binding 
sites to Au with the corresponding contact geometries called S-HH, S-BB, and S-TTl, respectively. For NH2 and CN we 
consider binding to single Au atoms in two different top positions with the contacts named NH2-TTI, NH2-TT2 and CN-TTf , 
CN-TT2. 
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35.1 0.002 


12.9 1.140 


17.9 0.907 


32.9 0.343 


33.2 


0.150 


33.1 


-2.389 


35.3 


-1.566 
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0.3 0.064 


0.1 1.313 


0.4 1.280 


0.2 0.436 


0.7 


0.201 


0.7 


-2.281 


0.1 


-1.423 


M2 


19.8 0.048 


19.3 1.208 


17.2 1.266 


20.4 0.429 


19.1 


0.191 


19.5 


-2.252 


20.2 


-1.404 


M3 


42.3 0.032 


42.1 1.127 


42.0 1.200 


46.3 0.382 


46.3 


0.189 


45.1 


-1.878 


46.4 


-1.158 


M4 


60.7 0.034 


53.0 1.006 


60.8 0.981 


61.6 0.400 


58.4 


0.158 


59.2 


-1.657 


59.9 


-0.938 


M7 


89.6 0.019 


84.0 0.298 


83.4 0.981 


87.4 0.191 


87.3 


0.091 


89.6 


-1.173 


89.9 


-0.632 



Table I: Torsion angle tp in units of degrees and the thermopower Q at T = lOK in units of ^V/K for all junction geometries. 



In the following the thermopower is calculated, if not 
otherwise indicated, for a low temperature of T = 10 K, 
where our theory is expected to apply best. For the 
DFT-based results presented below we determine Q by 
means of Eq. ([T|, i.e. by taking into account the thermal 
broadening of the electrodes. For the molecular junc- 
tions studied here, the differences to the values obtained 
via Eq. ^ often turn out to be small even at room 
temperature (T — 300 K). Hence thermopower values 
for higher T can be estimated using the values at 10 K 
through Q{T) w (T/IOK) x Q(IOK). Since we are not 
primarily interested in the temperature dependence of Q 
in this work, we suppress from here on the temperature 
argument. 



III. RESULTS AND DISCUSSION 

A. Models for the angle-dependent thermopower 

In Ref . we argued that the thermopower should de- 
pend on the torsion angle as 



Our argument was based on the observation that for a tt- 
orbital TBM in the off-resonant transport situation, the 
transmission of the bip henyl molecule can be expanded 
in powers of cos^ (f a<^12i2lIM\ 

T^{E) = a2iE) cos^ ip + a4{E) COS^ (fi + 0{cos^(p). (4) 



Conductance measurementsi^HiSland corresponding DFT 
calculations,'^ which both determine the transmission at 
the Fermi energy, show that 012 is the dominant term. 
The leading term in the tp dependence of Q is obtained 
from Eq. (|2| by taking into account the energy depen- 
dence of the expansion coefficients aj (E) and considering 
the terms up to j — 4. Then, we obtain Eq. ^ with 



a2{E) 



a2{E)dEai{E) - ai{E)dEa2{E) 



a2{Ef 



(5) 



.(6) 



E — EfT 



K, a + 6 cos p. 



(3) 



While this model uses minimal information about the 
biphenyl molecular junction, a disadvantage is that the 
magnitude and energy dependence of the coefficients a2 
and Oii are a priori unknown. 
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An alternative strategy is to use the 2LM of Ref. [T5] 
This minimal model explains the cos^ ^ law of the con- 
ductance by considering the pair of hybridizing frontier 
orbital resonances of the phenyl rings which are closest 
to Ep. Within this model, the transmission is given byiSi 



Yt cos Lp 



[E - isiv) - iT/2){E - - iT/2) 



(7) 



with a(iy9) = eo i tcoinp. Here, Eq is the relevant fron- 
tier molecular orbital energy of the individual phenyl 
ring. For the biphenyl molecule, it can be determined 
as the highest occupied molecular orbital (HOMO) or 
lowest unoccupied molecular orbital (LUMO) energy for 
vanishing inter-ring coupling t cos at (ys = 90° . The 
angle-dependent inter-ring coupling leads to a splitting 
of the pair of degenerate levels Eq at energies is^ai,^) — 
Eq i i cos Lp with symmetric and antisymmetric wave- 
functions, respectively. In addition, we have made the 
wide-band approximation with a symmetric and energy- 
independent coupling r to the left and right phenyl 
rings. The 2LM is hence characterized by the param- 
eters eo, r. 



We set e = — Ep, x = tcosif/^e^+T'^/'i, and 
assume |a;| <C 1. Performing a Taylor expansion in i, we 
obtain Eq. (Isl with 



4e 



a = —q- 



f2/4' 
4Pe(£2_3p2/4 

(e2+f2/4)^ 



(8) 
(9) 



These expressions predict that the sign of a,b is deter- 
mined by £. Thus, when e changes sign, a, b change 
sign at the same time. In the typical ofF-resonant trans- 
port situation |e| ^ |<| ,r, the sign of a w — (?4/e and 
b w —qAt^/e^ is identical. However, b may be of a dif- 
ferent sign than a in a more on-resonant case when the 
broadening F is of a similar size as e, i .e . , when — SF^ /4 
changes sign. 



B. Thermopower based on density functional 
theory 

For each of the biphenyl molecules in Fig. [T] we study 
the three different anchoring groups X = S, NH2, CN and 
select a total of seven contact geometries, as displayed 
in Fig. [2j F or S a nchors we choose three representative 
binding sites jlSSlMl where S binds covalently either to three 
Au atoms in the hollow position (S-HH), to two of them 
in the bridge position (S-BB) , or to a single one in the top 
position (S-TTl). NH2- and CN-terminated molecules 
bind selectively to a single Au electrode atom at each 
side via the nitrogen lone pair.^^''* Thus, we consider 




10 20 30 40 50 60 70 80 90 
(p (deg) 

Figure 3: (Color online) (a) Evolution of Q with increasing 
ip for all contact geometries. The symbols represent the ther- 
mopower values calculated with DFT at T = 10 K, and the 
lines are obtained by fitting Eq. Q to M1-M4 for each type 
of junction, (b) Zoom in on the Q values for S-HH, NH2- 
TTl, and NH2-TT2. (c) Schematic of the studied biphenyl 
derivatives and definition of the torsion angle p. 



two different top sites for NH2 (NH2-TTI, NH2-TT2) 
and CN (CN-TTl, CN-TT2), respectively. 

In Table [l] we summarize the torsion angle ip, which 
is defined as the dihedral angle between the two phenyl 
rings (see Fig.[3|, and the thermopower for all 42 molec- 
ular junctions studied. The data is presented graph- 
ically in Fig. |3] by plotting Q as a function of Lp for 
each of the seven types of junctions in Fig. [2] We no- 
tice that the sign of the thermopower is determined by 
the anchoring group. For the electron-donating S and 
NH2 linkers^Sl the energy of the 7r-electron system of 
the molecules is increased compared to the hydrogen- 
terminated case (X = H in Fig. [l|. The HOMO en- 
ergy is therefore close to Ep, as visible also from the 
transmission curves in Fig. |4j The hole conduction 
through the HOMO yields Q > 0, in agreeme nt w ith 
previous experimentai^ and theoretical results .'^^'i^ In 
contrast to this, for the electron-withdrawing CN an- 
chorin g group^ we have electron transport through the 
LUMoSEaHEll (see also Fig.|4|, and consequently Q < 0. 

Considering the absolute values of the thermopower. 
Fig. |3] shows that Q can differ markedly for the types 
of contact geometries. Given the off-resonant transport 
situation suggested by the transmission curves in Fig. [4] 
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and using Eqs. ([8| and (|9|, we can understand the results 
by changes in the level alignrnent e. As we will discuss 



in more detail below in Sec. IIIC level broadenings T 
and couplings t play no important role in that respect. 
The level alignment is determined by the charge transfer 
between the molecule and the electrodes, which is sen- 
sitive to the binding site of the anchoring group at the 
molecule-metal interface. For the thiolated molecules we 
find that the thermopower for S-BB and S-TTl is com- 
parable, but the values are significantly larger than those 
for S-HH. This behavior is related to our recent findings 
for the conductance of the thiolated molecules, where 
top and bridge geometries yield similar but much larger 
conductances than those with hollow sites.'^Sl Both obser- 
vations are due to a HOMO level which is more distant 
from Ep for S-HH as compared to S-BB and S-TTl. We 
explain this by the leakage of electrons from the molecule, 
including the S atoms, to the Au electrodes, when going 
from the S-TTl over the S-BB to the S-HH geometry. 
For the amines NH2-TTI gives a larger thermopower 
than NH2-TT2. We have checked that this is a result of 
the larger negative charge on the molecule when bonded 
in NH2-TTI position as compared to NH2-TT2, which 
moves the HOMO closer to Ep. With respect to the 
thiols we see that both NH2-linked geometries give rise 
to a thermopower well below those of S-BB and S-TTl 
but still larger than for S-HH. The CN- linked molecules 
show the largest \Q\. The more positive charge on the 
molecules in CN-TTl as compared to CN-TT2 leads to 
their smaller, i.e., more negative Q. 

Regarding MO with X = S, NH2 we can compare to ex- 
perimental and theoretical results for Q in the literature. 
For biphenyl-diamine a thermopower of Qmq^ — 
4.9±1.9 fiY/K was found at T = 300 K.lslwhich compares 

NH -TTl 

reasonably well to our calculated values of Q^o^' ~ 
10.52 /xV/K and Q^^~™ = 4.6/iV/K for the same T. 
Furthermore, recent calculations within a DFT approach 
with an approximate self-interaction correction for com- 
parable geometries showed similar results to ours.^^ For 
biphenyl-dithiol the comparison is complicated by the 
fact that our calculated values vary by two orders of mag- 
nitude for the different geometries, i.e., Qm™ — ^.11 
//V/K, Qfi^B ^ 39J4 ^v/K, Qf^JT ^ 28.08 ^iV/K at 
T = 300 K. They scatter indeed around the experimen- 
tal result of Q^^^P = 12.9 ± 2.2 /iV/K.STo our knowl- 
edge, the thermopower of cyano-terminated biphenyls 
has not yet been reported. A trend by the DFT calcula- 
tions to overestimate the thermopower can nevertheless 
be recognized.!^ It is expected from the typical overesti- 
mation of experimental conductance values, attributed 
mostly to the interpretation of Kohn -Sham eigenvalues 
as approximate quasi-particle energies. fSSEIl According to 
Eqs. ([8| and ^ an underestimation of |e| leads to an 
overestimation of \Q\. However, finite temperature ef- 
fects due to vibrations, not accounted for in our calcu- 
lations, may also play a role in the room-temperature 
experiments. 

The transport through the well-conjugated molecules 





a (mV/K) 


b (^V/K) 


r (%) 


S-HH 


0.02 


0.04 


31 


S-BB 


0.86 


0.43 


13 


S-TTl 


0.93 


0.37 


13 


NH2-TTI 


0.37 


0.06 


4 


NH2-TT2 


0.15 


0.05 


12 


CN-TTl 


-1.44 


-0.89 


16 


CN-TT2 


-0.81 


-0.65 


20 



Table H: Parameters a and b of Eq. ([3| used in Fig. |3]to fit 
the DFT results of M1-M4 for each type of junction, and the 
relative change r of Q between Ml and M4. 



M0-M4 is dominated by the tt electrons, and we have 
shown in Refs. I18I19I281 that for these molecules the 
conductance arises from one transmission eigenchannel 
of TT character .liSllll Hence we would expect their ther- 
mopower to follow Eq. ([3|. Despite the variations of Q 
with anchoring groups and binding positions, we find a 
weak cos^iy9-like decrease of the absolute values for Ml- 
M4 for all types of geometries. MO, however, deviates 
from this trend. Although the electron-donating effect of 
the alkyl chains is expected to be small, it increases Q 
for M1-M4 as compared to MO. To clarify this we cal- 
culated by means of electrostatic potential fitting and a 
Lowdin population analysis the charge transferred from 
the alkyl side chains to the two phenyl rings for the 
hydrogen-terminated (X = H in Fig. [ij, isolated gas- 
phase molecules. Both methods yield an overall negative 
charge on the phenyl rings which is practically indepen- 
dent of the alkyl chain length. Therefore the substituent- 
related energy shift of frontier orbital levels is similar for 
M1-M4, and the a + 6cos^ (p dependence is observed. 

Focusing on the thermopower of M1-M4, we extract a 
and b by fitting their Q with Eq. (|3|. The precise values 
are given in Table [TTj and the corresponding fits are shown 
as continuous lines in Fig.[3j Additionally, we list in Table 

pjtheratior = |Qm1-'3m4I/('3m1+<5M4)' ^^^^tify^S 
the maximal decrease of Q in that subset of molecules. 
We find it to vary between 4% and 31%. In detail, we 
observe the largest relative change for S-HH followed by 
CN-TT2. CN-TTl, S-BB, S-TTl, andNH2-TT2 all show 
a similar r, while it is smallest for NH2-TTI. 

For M7, Eq. ([3| is not expected to hold, because 
the transport at ~ 90° is not 7r-like but proceeds 
through transmission eigenchannels of n-a character.'^^'^ll 
Furthermore, M7 shows the largest substituent-related 
shifting effect on the biphenyl backbone in our fam- 
ily of molecules due to the electron-d onating nature of 
the four attached methyl side-groups. I^ESIiQl Jtg ther- 
mopower hence arises from a detailed interplay between 
the substituent-related shifting and the large torsion an- 
gle, as explained in Ref. 22 We find that the absolute 
values of Q for M7 are generally lower than predicted by 
the fits with Eq. ^. Only for S-TTl and NH2-TT2 the 
thermopower seems to follow the a + b cos^ (p dependence 
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but this is likely coincidental. 



C. Transport analysis using the vr-orbital 
tight-binding model and the two-level model 

In order to better understand the differences in the 
thermopower for the various anchoring groups and bind- 
ing positions, we need to examine the parameters e,t,T 
of the 2LM which determine the thermopower according 
to Eqs. ([8| and (|9|. We note that the dominant, angle- 
independent term a is a function of e, T only. Thus, to 
discuss main anchor-group- and binding-site-related vari- 
ations of Q for the seven different junction types of Fig.|2] 
it is sufficient to concentrate on these two parameters. 
The If dependence of Q, however, results from the in- 
terference of the hybridizing pair of phenyl-ring frontier 
orbital levels, and b hence depends also on t. 

We obtain the parameters of the 2LM from the TBM 
introduced in Ref. :32] The TBM is sketched in Fig. ga). 
Similar to the 2LM, the Hiickel-hke TBM is character- 
ized by three parameters which are the on-site energy 
£o of each carbon atom, the nearest-neighbor hopping 
t between atoms on each of the phenyl rings, and the 
electrode-related broadening F. The inter-ring hopping 
is given as t' — t cos ip. Using the wide-band approxima- 
tion, we assume all components of the lead self-energy 
matrices to vanish except for = (S/j)^^ = — iF/2, 

with a and w indicating the terminal carbon atoms of the 
biphenyl molecule as shown in Fig. |4|^a). 

The parameters eo,t,T of the TBM are extracted by 
fitting t{E) curves calculated with DFT. We focus on 
the molecules M1-M4 and set ip to the torsion angle re- 
alized in the specific junction geometry (see Table ^. 
Concentrating particularly on the HOMO-LUMO gap 
and frontier orbital peaks, we find that the fitted TBM 
generally reproduces well the transmission in that range 
and that the parameters extracted for M1-M4 are very 
similar in each of the seven types of junctions. Fi- 
nally, the parameters of the 2LM are derived from those 
of the TBM. e and i are obtained by evaluating ap- 
propriate eigenvalues of the angle-dependent Hiickel-like 
Hamiltonian of the TBM. For F we identify imaginary 
parts of complex eigenvalues of the non-Hermitian ma- 
trices {H + T,l + for the TBM and the 2LM, 
respectively.'^ Here, Hjk and represent 
the matrix elements of the Hamiltonian and of the elec- 
trode self-energies in the corresponding model. All the 
parameters determined in this way are listed in Table [ITT] 
For M2, transmission curves calculated with the DFT, 
the TBM, and the 2LM are shown in Fig. gb) for each 
of the three anchoring groups. 

Using the parameters of Table |ni| we compare in 
Fig. gc) Q as a function of (f for the TBM and 2LM 
fits with the DFT results. We find that the TBM agrees 
well with the DFT-based values for the illustrated junc- 
tion geometries S-TTl, NH2-TTI, and CN-TT2. The 
2LM, instead, overestimates |Q| somewhat. Considering 





eo 


t 


r 


e 


t 


r 


S-HH 


-4.40 


-2.30 


0.70 


-1.70 


-0.68 


0.22 


S-BB 


-4.02 


-1.95 


1.10 


-0.97 


-0.58 


0.35 


S-TTl 


-4.00 


-1.90 


0.96 


-0.90 


-0.56 


0.31 


NH2-TTI 


-4.30 


-2.29 


0.60 


-1.59 


-0.68 


0.19 


NH2-TT2 


-4.40 


-2.32 


0.66 


-1.72 


-0.69 


0.21 


CN-TTl 


-6.10 


-2.00 


0.14 


0.90 


-0.59 


0.04 


CN-TT2 


-6.05 


-1.99 


0.15 


0.94 


-0.59 


0.05 



Table III: Parameters of the TBM eo, t, ^ obtained by fitting 
the DFT-based transmission curves for M1-M4 for each type 
of junction. The parameters e, t, V of the 2LM are derived 
from those of the TBM as described in the text. All values 
are given in units of eV. 



Eq. ^ and the transmission curves in Fig. gb), we at- 
tribute this to an underestimation of t{Ep) and an over- 
estimation of \dET{EF)\. All results exhibit a consistent 
weak dependence of Q on ip. 

The data in Table |fll| shows that transport through 
the biphenyl molecules is off-resonant with the re lation 



III A 



F ^ |e| being well fulfilled. As argued in Sec 
a, b should thus take the same sign and change it to- 
gether with Q when the transport for S- and NH2-linked 
molecules changes from HOMO- to LUMO-dominated 
for CN anchors. This is consistent with our findings in 
Figs, [sjandgc), and explains the decrease of |Q| with 
increasing ipW^ 

Coming back to our discussion of the differences of the 
thermopower for the various anchoring groups and bind- 
ing positions in Fig. |3] we observe that e is around 0.6 
to 0.8 eV closer to Ep for S-BB and S-TTl as compared 
to S-HH, NH2-TTI, and NH2-TT2, which explains their 
larger Q. For the CN-terminated molecules |e| is com- 
parable to those for S-BB and S-TTl. SHghtly larger 
values of IQI for CN result from the very small broaden- 
ings F. Furthermore, for both NH2 and CN, i and F are 
essentially independent of the binding position, and the 
difference in Q between TTl and TT2 hence stems from 
the changes in the alignment of the HOMO and LUMO 
levels. 



IV. CONCLUSIONS 

We have analyzed theoretically the thermopower of 
single-molecule junctions consisting of biphenyl deriva- 
tives contacted to gold electrodes. Or DFT-based study 
with the three anchors S, NH2, and CN shows a positive 
thermopower for S or NH2 and a negative one for CN. 
For the junction geometries considered, different binding 
sites did not affect the sign of Q but led to variations 
of absolute value. For thiolated molecules in bridge and 
top binding sites Q can be up to two orders of magnitude 
larger than for molecules bonded in hollow position, while 
the variations for the two considered top binding sites 



7 




* 0.5 
I " 

-1 

-1.5 
-2 



— 1 


1— ' — 1 — ' — I 

• T ■ 


— 1 — 1 — ^ 


1 P 1— 1- 

— ^ ^ ^ " 


—1 — ' — r- 


1 — 1 — ' — 1 — ; 




♦ 








* 




























S-TTI 


• DFT 


— TBM 


--2LM 






NH,-TTI 


» DFT 


— TBM 


- - 2LM 






CN-TT2 


DFT 


TItM 




















1 1 1 , 1 


1 1 L. 


1 1 1 i_ 


_i 1^ i_ 





10 20 30 40 50 60 70 80 90 

Figure 4: (Color online) (a) Schematic of the TBM used to fit 
DFT-based transmission curves, (b) Transmission of M2 as 
a function of energy calculated with DFT, and the fits using 
the TBM and the 2LM. (c) Q as a function of i/J, comparing 
values obtained with the TBM and the 2LM to the DFT- 
based results. In panels (b) and (c) S-TTl, NH2-TTI, and 
CN-TT2 junction geometries were selected. 



were around a factor of two for NH2 and CN anchors. 
We have explained these observations by the changes in 
the level alignment of current-carrying frontier molecular 
orbitals. They are caused by the binding-site-dependent 
charge transfer at the metal-molecule interface. 

The main purpose of this work was the study of the 
dependence of the thermopower on conjugation for an ex- 
perimentally relevant system. In our set of six biphenyl 
derivatives, the conjugation was controlled by the torsion 
angle </? between the phenyl ring planes, and it was varied 
stepwise between and 90° by means of alkyl side chains 
attached to the molecules. Despite the sensitivity of the 
thermopower to the precise geometry at the molecule- 
metal interface, we observed for all investigated types of 
junction configurations a decrease in \Q\ with increas- 
ing (/?, following a characteristic a -I- &cos^ ip law. We ex- 
plained this behavior in terms of a two-level model, which 
considers the pair of hybridizing frontier orbitals on the 
phenyl rings. Predictions by this model of a simultane- 
ous change in sign of Q,a,b for a change from HOMO- 
to LUMO-dominated transport in the ofF-resonant situ- 
ation are consistent with our DFT results. Overall, the 
influence of conjugation on the thermopower is much less 
pronounced than on the conductance. 

We propose to measure the a + b cos^ (p dependence of 
the thermopower for the set of biphenyl molecules studied 
here. Using alkyl chains of different lengths, parasitic 
substituent-related shifts in Q, superimposed on the weak 
a + &cos^ ip dependence, are largely avoided. Depending 
on binding site and employed anchoring group, relative 
variations of Q of around 15% are expected between Ml 
and M4. Since frontier molecular orbital energies are 
likely positioned closer to Ep in our calculations than 
in the experiment, the relative changes of Q with ip are 
expected to be somewhat smaller than in our theoretical 
predictions. Nevertheless, we suggest that the variations 
of the thermopower with torsion angle are experimentally 
detectable. 
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